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Abstract. Classical theories of turbulence do not describe accurately inertial range scaling laws 
\Q . in turbulent convection and notably fail to model the shape of the turbulent spectrum of solar 
photospheric convection. To understand these discrepancies, a detailed study of scale-by-scale 
budgets in turbulent Rayleigh-Benard convection is presented, with particular emphasis placed 
on anisotropy and inhomogeneity. A generalized Kolmogorov equation applying to convection is 
^ • derived and its various terms are computed using numerical simulations of turbulent Boussinesq 
convection. The analysis of the isotropic part of the equation shows that the third-order velocity 
structure function is significantly affected by buoyancy forcing and large-scale inhomogeneities. 
Anisotropic contributions to this equation are also shown to be comparable to their isotropic 
counterpart at moderate to large scales. Implications of these results for convection in the solar 
photosphere, mesogranulation and supergranulation are discussed. 
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<n : 

OO ■ 1. Introduction 

Fundamental information about the physical processes at work in turbulent flows can 
\q | be obtained by studying their energy spectrum and the associated spectral scaling laws. 
C^) ■ Various studies of the solar surface have focused on understanding the spectrum of tur- 
bulent convection in the photosphere (Zahn 1987; Espagnet et al. 1993; Hathaway et al. 
2000). In spite of these efforts, a globally satisfying theoretical picture of the distribution 
of energy at various scales (like granulation, mesogranulation, supergranulation) is still 
lacking (Rincon, Lignieres & Rieutord 2005): solar convection does not seem to follow 
the simple predictions of linear stability theories or those of classical theories of isotropic 
C$ • and homogeneous turbulence. 

There are several possible explanations for this disagreement between classical theories 
• ■ *h . and observations (see for instance Petrovay 2001). On the one hand, solar convection is 
not laminar, therefore linear theories will fail unless mean-field effects (like turbulent 
viscosity) can be parametrized accurately. On the other hand, turbulence in the solar 
photosphere is neither isotropic (because it is driven by the unstable stratification in the 
vertical direction) nor homogeneous (because the observed flow occurs in the vicinity of 
the t = 1 optical depth surface, which acts as a natural boundary). Buoyancy effects 
are also very likely present at most observable scales, which are thus not strictly part 
of an "inertial" range in the Kolmogorov sense. The focus of this essay is on these 
last three points. Results regarding inertial range scalings in numerical simulations of 
turbulent convection, which have been obtained using modern tools of turbulence theory, 
like the SO (3) decomposition of structure functions to quantify anisotropic effects and a 
description of statistical inhomogeneities, are used to illustrate how convection spectra 
are affected simultaneously by forcing, anisotropic and inhomogeneous effects. 
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In §[21 a quick summary of the current status of theories of turbulent convection is 
given. §[3] is devoted to the analysis of turbulent statistics and scaling laws in a numerical 
simulation of turbulent Boussinesq convection. A discussion of these results in the context 
of the solar photosphere is presented in §[?]and a short conclusion follows. 



2. Theories of turbulent convection 

The aim of this section is to present a quick summary of currently available theories 
to describe turbulent convection. Emphasis is placed on a description in physical space 
using structure and correlation functions rather than on turbulent spectra, because scal- 
ing predictions result from exact laws (Kolmogorov and Yaglom equations, see Monin & 
Yaglom 1975) which are valid in physical space only. Intermittency effects are not con- 
sidered. Moments of increments of various quantities such as velocity u or temperature 9 
fluctuations between two points separated by a correlation vector r are used, and incre- 
ments of any quantity / are denoted by Sf. Simple forms for these exact laws are usually 
obtained assuming homogeneity (turbulence looks statistically the same everywhere) and 
isotropy (turbulence looks statistically the same for any orientation of r). 

2.1. Kolmogorov theory 

Kolmogorov's 1941 theory (hereinafter K41) predicts 5u ~ r 1 / 3 and 59 ~ r 1 / 3 . The first 
prediction can be obtained using dimensional analysis or one of the only exact results 
known for turbulence, namely the Kolmogorov equation, also called the 4/5 law: 

((6u r f) = ~ (e) r + 6^ <(^ r ) 2 ) , (2.1) 

where v is the kinematic viscosity and (e) = v (ipiUj) 2 ). This law states that the third- 
order structure function for longitudinal velocity increments ((fe r ) 3 ) is proportional to 
r in the inertial range of turbulence, where both viscosity and forcing are negligible. The 
prediction for temperature increments (at scales larger than the thermal diffusion scale) 
results from the extra assumption that temperature behaves as a passive scalar, which 
is strictly speaking only valid for a neutrally stratified flow. The corresponding spectra 
are both proportional to fc~ 5 / 3 in the framework of this theory. 

2.2. Bolgiano-Obukhov theory 

Atmospheric science research has led to the development of Bolgiano-Obukhov (Bolgiano 
1959 ; Obukhov 1959, hereinafter B059) theory for stably stratified turbulence. This 
theory is equally valid for an unstably stratified medium (L'vov 1991). Unlike K41, it 
does not require that the forcing acts on large scales only and relies on a dominant 
balance in a "generalized" 4/5 law between ^((5w r ) 3 ) and a forcing/damping term due 
to unstable/stable stratification (Yakhot 1992), whereas in the classical 4/5 law, the 
balance is between ^(5u r ) 3 ') and the 4/5 term. Using the equivalent of the 4/5 law for 
temperature, called Yaglom equation or 4/3 law, one can infer the following scalings 
using dimensional analysis: Su ~ r 3 / 5 and 59 ~ r 1 / 5 , leading respectively to fc~ n / 5 and 
fc~ 7 / 5 spectra for the velocity and temperature fluctuations. The crossover scale between 
K41 and B059 is called the Bolgiano length (Chilla et al. 1993; Benzi et al. 1998) 

Nv}/ 2 d 

LB = (RaPry/* ' (2 ' 2) 

where Nu is the Nusselt number, d is the typical scale-height of the layer, Ra is the 
Rayleigh number and Pr is the Prandtl number. At this scale, the 4/5 term equals the 
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forcing/damping term in the Kolmogorov equation. B059 should hold for r > Lb, while 
K41 predictions should be observed for r < Lb, leading to a priori observable kinks in 
both velocity and temperature spectra. Note that equation (|2.2|) is only a dimensional 
estimate which, from previous studies (Calzavarini et al. 2002, Rincon 2006), appears to 
underestimate the actual by a factor ~ 10. 

2.3. Dealing with anisotropy and inhomogeneity 

As mentioned earlier, important differences between theories and observations or exper- 
imental convection (e. g. Verzicco & Camussi 2003) may be related to anisotropic or in- 
homegeneous effects (Hill 1997). A brief description of modern techniques used to quantify 



these effects is given here. Many details on these issues are given by Biferale & Procaccia (2005) 



Danaila et al. (2001) and their applications to turbulent convection can be found in the 



work of Biferale et al. (2003) and Rincon (2006) The essential point is to project the full 



(x,r) dependent (as opposed to a |r| only dependence, with x being the measurement 
position) Kolmogorov and Yaglom equations onto a spherical harmonics basis. Projec- 
tions on the t = degree lead to the usual isotropic Yaglom and Kolmogorov equations. 
Inhomogeneous effects are included by keeping all terms involving variations of statistical 
quantities with respect to the location at which they are computed (e. g. at the centre). 
The isotropic part of the generalized Kolmogorov equation, used in §[3j reads 

- (U R )° a + 2 -^Jy 2 {Su z S9)° Q dy + 2v^-({8u l f)l + (NH)° Q = ^Jy 2 [(e) + (e')}° dy ,(2.3) 

where the [J notation indicates projection onto the Yq spherical harmonic (equivalently 
average over a r-sphere), (E/r)o = ((8ui) 2 5u r ) is a generalized third-order velocity 
structure function, (e') = (e) (x + r), the associated (e) + {e') term is a generalized form 
of the 4/5 term in equation (|2.ip . the ag term is the contribution of buoyancy forcing in 
the generalized Kolmogorov equation, and (NH) stands for all inhomogeneous terms. All 
these statistical quantities depend on the vertical (gravity) coordinate. This equation is 
exact: only statistical stationarity has to be assumed to obtain it. Evaluating its various 
terms to determine the dominant balances and thus potential inertial range scalings for 
(Ur} and (0/{)q is the aim of the next section. 



3. Numerical results 

Three-dimensional direct numerical simulations of turbulent Boussinesq convection at 
Ra = 10 6 and Pr = 1 with aspect ratio A = 5 have been performed in order to compute 
the scale-by-scale budgets of equation (|2.3p in the middle of the convection box. 

3.1. Isotropic component 
Figure [1] shows that the buoyancy term in Kolmogorov equation is smaller than the 4/5 
term, so that the conditions for B059 are not met. The two curves seem to intersect at 
scales comparable to the box depth (r/rj = 60, where r\ is the Kolmogorov dissipation 
scale) rather than 0.1 d, predicted by equation (|2.2[) for this simulation. However, K41 is 
not verified either. The reasons are that inhomogeneous and buoyancy terms interfere in 
the scale- by- scale budget to prevent the K41 balance between (Ur) and the 4/5 term. A 
detailed investigation (Rincon 2006) shows that Ra = 10 9 is a mimimum requirement for 
buoyancy to be negligible for some r < Lb and to observe one K41 decade. The reason 
why scaling laws are not observed on velocity spectra or structure functions is therefore 
very different from the usual argument that viscous dissipation (plotted in figure [1]) is too 
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large. In fact, scaling laws are observed in this simulation for mixed velocity-temperature 
third-order structure functions appearing in Yaglom equation. 




1 10 

Figure 1. Scale-by-scale budget for the isotropic part of Kolmogorov equation: log-log plots 

of - (Ur)° (- -), 2ag/r 2 £ y 2 (6u*S0)° o dy (- - -), 2vd r <(^) 2 )o (" ■ (NH)° ( ). The 

sum of the l.h.s. terms of equation (|2.3[) ( ) matches the r.h.s. "4/5 term" (■ ■ ■). 

3.2. Anisotropy and projection effects 

The other effect which may prevent observations of scaling laws is the anisotropic driving 
mechanism, which generates I ^ components of structure functions comparable to the 
isotropic component at moderate to large scales (figure [2] depicts the spherical harmonics 
spectrum of {Ur) for even £ and various r) . This is important because "reduced" structure 
functions based on correlation vectors lying in horizontal planes are linear combinations 
of these anisotropic components, which may follow different r-power laws according to 
modern turbulent theories (Biferale & Procaccia 2005): 

(U R ) (r, tt/2) = J2 (Ur)1 ■ (3.1) 

e 

Thus, unless turbulence is really isotropic, these reduced quantities behave differently 
from isotropic structure functions, and may not exhibit any scaling behaviour. The local 
scaling exponents of (Ur) (r, tt/2) and of the r.h.s. of equation (|3 . If) reconstructed up to 
£ = 6 arc plotted in figure [2] using the associated ^-spectrum. It is straightforward to 
see that the £ — 2 and £ = 4 degrees strongly affect the local slope of the "reduced" 
third-order structure function (slope 1 corresponds to the isotropic K41 scaling law). 

4. Understanding the spectrum of photospheric convection 

This section discusses the theoretical arguments exposed in §[2] in the context of the 
solar photosphere, in the light of the numerical results presented in §[3j 

4.1. Inertial and injection ranges for the photosphere 

Recent studies indicate that Nu ~ i?a 1 / 2 Pr 1 / 2 holds for Ra > 10 11 (Grossmann & Lohse 
2000; Calzavarini et al. 2005), so that Lb as given by equation (|2.2|) is independent of 
any nondimensional parameter in the hard turbulence regime and is of order a few d for 
large Ra (taking into account the prefactor). Close to r = 1, identifying d with the local 
density scale height H p leads to an actual Lb of a few thousands kilometers. Thus, surface 
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Figure 2. Left: spherical harmonics spectrum of (Ur) at the centre of the box for r jr\ = 7.9 

( ), r/ri = 15.3 (■ ■ ■) and r/n = 22.6 ( ). Right: local log-slope of (U r ) (r,n/2) ( ) and 

using spherical harmonics up to £ = (+), £ = 2 (x), i — 4 (o), £ = 6 (A) (see equation (13. ip ). 



turbulence from granular to supergranular scales is strongly affected by buoyancy forcing, 



and such scales are part of the injection range, as suggested earlier by Zahn (1987) and 



Petrovay (2001) This is supported by spectral-space budgets computed for large aspect 
ratio simulations by Rincon et al. (2005)| in which most energy is pumped predominantly 



at mesogranular scales and cascaded down to the dissipation scales. 

4.2. The effects of inhomogeneity and anisotropy on data analysis 

According to this estimate for Lb, K41 should be observed at subgranular scales, while 
B059 is expected at larger scales. In practice, there is no clear evidence for this. A reason- 
able explanation for this disagreement on large scales is that inhomogeneities dominate 
the scale- by- scale budget at scales larger than the typical scale-height and prevent the 
B059 balance (see Danaila et al. 2001 for large-scale inhomogeneous effects in a channelt 
flow) . As granular scales are close to Lb, there is no dominant balance in the Kolmogorov 
equation for this range of scales and therefore no definite inertial range scaling exponent. 
Finally, photospheric convection spectra presented so far in literature have been obtained 
using either line-of-sight measurements or by analysing local horizontal motions, i. e. 
along specific directions only. As shown previously, these quantities involve a mixture of 
different power laws and attempts to extract (isotropic) exponents from them may prove 
impossible. To disentangle anisotropic effects would notably require the knowledge of the 
vertical velocity field at various depths to measure correlations along the vertical. Only 
very high-resolution local helioseismology may provide this information in the future. 



5. Conclusion 

This study is only a first step towards a fully consistent description of the spectrum of 
photospheric convection, but shows that considering generalized Kolmogorov and Yaglom 
equations is crucial to understand anisotropic and inhomogeneous turbulence forced in 
a full range of scales. An interesting result is the independence of the estimated Lg on 
both Pr and Ra at Ra > 10 11 , which makes it directly applicable to the Sun. More work 
is needed to understand scaling laws in the low Pr regime typical of the Sun. 
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Discussion 

J. P. Zahn: Is your statement that energy is injected on mesogranular scales based on 
numerical simulations ? 

F. RlNCON : Not only. This is supported by theory (equation (|2,2[) with Nu ~ (RaPr) 1 ^ 2 ). 

F. H. Busse: What boundary conditions did you use in your Boussinesq simulations ? 

F. Rincon: Fixed temperature and stress-free boundaries are used for both plates. 

H. G. Ludwig: Can a stationary theory (like Kolmogorov's) describe the properties of 
a transient flow like photospheric convection? 

F. RlNCON: Yes. Statistical stationarity does not prevent time-dependence in the flow. 

H. G. Ludwig: Do you associate the observable granular cells with turbulent eddies ? 

F. Rincon: Granules are undoubtedly turbulent. But they are very much influenced by 
anisotropy, inhomogeneity and injection, unlike classical inertial eddies. 

R. F. Stein: On the Sun, large temperature fluctuations are at granular scales, which 
is where there are the large entropy fluctuations that give rise to the buoyancy driving. 

F. RlNCON: True. Granular scales are very much influenced by buoyancy at r = 1. 
Deeper and deeper, driving occurs on larger and larger scales, as Lb ~ a few H p shows. 



